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RESEARCH SUMMARY 


Equations were derived for predicting crown width of trees 
from diameter, height, crown length, and basal area per acre, 
and for predicting foliage weight of trees from diameter, height, 
crown length, age, relative diameter, and number of trees per 
acre. Coefficients were estimated for 11 conifer species in 
northern Idaho and western Montana. Embedding these equa- 
tions in the prognosis model for stand development will en- 
hance the prediction of vegetation characteristics needed for 
interfacing insect outbreak, wildlife habitat, and watershed 
models. 
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INTRODUCTION 


This paper reports on the development of multiple linear 
regression models predicting individual tree crown width and 
foliage weight for 11 conifer species in northern Idaho and 
western Montana. The equations have been applied in a spruce 
budworm version’ of the stand prognosis model (Stage 1973) 
being currently developed. The crown width and foliage bio- 
mass equations, when linked with, other relationships in the 
spruce budworm system, help predict insect damage and dis- 
persal of larvae during simulated outbreaks. The crown width 
function has also been used to predict percent stand crown 
cover in the prognosis model (Stage 1973) for linkage with 
proposed watershed’ and big game habitat prediction® sub- 
models. 

The relationship between crown width and tree diameter has 
been of interest to several researchers in developing crown 
competition indices (see for example Krajicek and others 1961 
or Vezina 1963). This relationship has been expressed as: 


crown width = a + b (d.b.h.) 


‘Crookston, N.L. The interface between the stand prognosis model and the 
western spruce budworm model — Version 1.2. Unpubl. rep. on file, Intermt. For 
and Range Exp. Stn., Forestry Sciences Laboratory, Moscow, Idaho, October 
1980, 30 p. 

Simons, D.B., R.M. Li, and T.J. Ward. 1980. Development of a generalized 
planning mocel for assessing effects of land use:changes on watershed re- 
sponses and aquatic habitat conditions. Unpubl. study plan, 14 p. Colo. State 
Univ. Dep. Civil Eng., Fort Collins. 

3Peek, J.M. 1970. Study plan: evaluations of elk habitat use patterns and 
relationships in the Gospel-Hump area of central Idaho, 10 p. Univ. Idaho, Wild! 
Resour., Dep. For., Wildl. and Range Sci., Moscow. 
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In some studies, for example Stiell (1966), stand density in- 
versely influenced tree crown width. Vezina (1964) used stand 
basal area to predict stand crown closure. Bonnor (1964), in 
studying the relationship between diameter at breast height 
(d.b.h.) and crown width on vertical photography, found no 
effect from stand density. The current study investigates the 
contribution of a number of stand and individual tree character- 
istics in predicting crown width by species. 

The second part of this study looks at the relationship be- 
tween individual tree foliage weight and measured tree and 
stand characteristics. Kittredge (1944) developed the following 
logarithmic expression to estimate leaf weight for several conif- 
er anc hardwood species: 


In (leaf weight) = b (In d.b.h.) — a 


He noted conflicting site and age effects, possibly confounded 
by geographical location. The effects of density and crown class 
as sources of variation were not separated out in his study. 

Brown (1978) approached the problem of estimating foliage 
weight indirectly by first developing total crown weight as a 
function of tree diameter, height, and crown length, and then 
estimating the proportion of the total crown in foliage in a sepa- 
rate equation. In addition, Brown presented more than one 
model form for each species depending on crown class and size 
class. The current study used Brown's data with the purpose of 
simplifying the direct estimation of foliage weight by investi- 
gating one or two model forms applicable to al! tree species 
considered. 


DATA 


The data studied are a subsample of trees used by Brown 
(1978) in estimating crown weights. Trees deleted from the 
original data set were those lacking measurements for crown 
width or foliage weight. Refer to Brown (1978) for a more 
complete description and listing of the data. 

Measurements from dominant and codominant trees were 
studied for the following 11 species (abbreviations are used in 
following tables and figures): 


DF Douglas-fir Pseudotsuga menziesii (Mirb.) 
Franco 

ES Engelmann spruce Picea engelmannii Parry 

GF Grand fir Abies grandis (Dougl.) Lindl. 

LP Lodgepole pine Pinus contorta Doug. 

PP Pondercsa pine Pinus ponderosa Laws. 

AF  Subalpine fir Abies lasiocarpa (Hook.) Nutt. 

WH_ Western hemlock Tsuga heterophylla (Raf.) Sarg. 

WL Western larch Larix occidentalis Nutt. 

WC_ Western redcedar Thuja plicata Donn 

WP Western white pine Pinus monticola Dougl. 

BP Whitebark pine Pinus albicaulis Engelm. 


In addition, trees of intermediate and suppressed crown classes 
were available only for western redcedar (WC-INT), ponderosa 
pine (PP-INT), grand fir (GF-INT), and Douglas-fir (DF-INT). 
Tables 1 and 2 show the distribution of data used in developing 
the models. 


PROCEDURE 


The models were developed by means of a regression 
screening method. First, the relationships between the two 
dependent variables of crown width and foliage weight, and 
stand and tree characteristics were studied by screening for the 
best combinations of the independent variables and their log- 
arithmic transformations. Variable screening and coefficient 
estimation were done in a multiple linear regression computer 
program, REX (Grosenbaugh 1967). The following relation- 
ships between the dependent and independent variables were 
examined: 


Y = {(0) D2,,H, CL\CR,A, TPA, BA DREL SSP. Si) 


InY = f(InD, (InD)?, InH, INCL, InCR, InA, InTPA, 
InBA, InDREL, SP, SI) 


Table 1.—Distribution of sample trees by species and crown class used in 


estimating crown width 


Number 
of trees 
<3.5 23.5 

Species' inches d.b.h. d.b.h. 

Inches 
BP 12 1 0 - 10.6 
WP 17 8 0 - 24.7 
WC 14 14 0 - 26.6 
WC-INT 11 11 0 - 10.6 
WL 15 9 0 - 21.8 
WH 16 12 0 - 21.4 
AF 15 10 0 - 12.7 
PP 15 20 0 - 34.0 
PP-INT 12 9 0-12.14 
iP 12 7 0 - 15.6 
GF 16 13 0 - 20.4 
GF-INT 13 fe) 0- 11.8 
ES 14 8 0 - 23.2 
DF 16 13 0 - 33.9 
DF-INT 12 10 0-11.2 


Total 210 160 


'See text page 2 for species list. 


Range of Measurements 


Crown Stand 
Height length basal area 
weceeenennnene Feet --------------- Feet?/acre 
3- 44 2- 36 1 - 181 
3 - 126 3- 88 1 - 299 
2-127 298 1 - 300 
4- 64 2- 42 45 - 335 
3-129 Shot 7/s) We 227A 
4-126 3-111 1 - 320 
3- 93 2- 70 1 - 260 
2-144 2- 89 1 - 325 
4- 70 1- 44 22 - 280 
2- 84 1- 60 Le ZA 
3-137 3 - 108 1 - 432 
4- 72 2- 60 39 - 296 
2-128 2-105 1 - 240 
4-143 4- 86 2 - 238 
4- 83 lier 477 36 - 360 


Table 2.—Distribution of sample trees by species and crown class used in estimating 
foliage weight 


SS eS eS eee 


Range of measurements 


Number 
of trees Number 
<35) S35 Crown of trees 
Species' inches d.b.h. d.b.h. Height length per acre Age 
Inches _ -------------- Feet ------------- Years 
BP 12 4 O- 7.4 3- 38 7, 461 - 6324 6 - 182 
WP alg/ 2 O= 7:4 3- 66 3-46 300 - 4500 4- 55 
WC 14 11 0 - 26.6 2-127 2 - 98 110 - 7800 3 - 237 
WC-INT 11 11 0- 10.6 4- 64 2-42 169-3611 16-162 
WL 15 3 O- 6.6 3- 50 3-29 600 - 5100 S155 
WH 16 4 O- 7.0 “A= 38 3 - 38 848 - 7200 AS Se 
AF 15 10 0- 12.7 3- 93 2-70 87 - 15121 4-201 
PP 15 20 0-340 2-144 2-89 5 - 7200 6 - 217 
PP-INT 12 9 0- 12.1 4- 70 1-44 142 - 3344 11-140 
LP 12 6) O-gowe 2.-" 413 eons 2 300 - 6000 Se ail 
GF 16 8 O-15:6 5 3= 79 3-70 300 - 7079 4-124 
GF-INT 13 9 O-J1e8 “Ai. 72 2-60 359-5804 15-117 
ES 14 s' O- 9.0 2- 57, 2-44 600 - 12730 7 - 153 
DF 16 10 0-339 4-143 4-86 21 - 7200 6 - 262 
DF-INT ie 10 0- 11.2 4- 83 1-47 261-4500 12-145 
Total 210 114 
‘See text page 2 for species list. 
where: 
RESULTS 


D = diameter breast height, inches (cm) 
H = height, ft (m) 

CL = crown length, ft (m) 

CR = crown ratio 
A = tree age, years 

TPA = number of trees per acre (number per hectare) 
BA = stand basal area, ft?/ac (m*/ha) 
DREL = relative diameter (tree d.b.h/quadratic mean 
stand diameter) 


1 if tree is of species i, 

0 otherwise 

1 if tree is in site class j, 
Sl; = 


0 otherwise 


Y = predicted crown width (CW), ft (m) or 
foliage weight (WT), Ibs (kg) 


Variable selection began with the model chosen in the REX 
screening process as having the lowest mean square error 
var(Y). A general model was chosen when coefficients that 
were not significantly different from zero when tested by the 
T-statistic at the 0.05 level, were deleted from the parameter 
set, and all remaining coefficients were significant when the 
model was refit. 

Refinements on the general model were investigated by 
allowing each parameter to vary by species. Species x para- 
meter interactions that significantly improved the regression 
sum of squares were then included in the final model form. 

In developing the equations, the data were divided into two 
sets consisting of those trees 3.5 inches (8.9 cm) d.b.h. and 
larger, and trees less than 3.5 inches (8.9 cm). Each set was 
then modeled independently. The break in diameter was 
chosen to correspond to prognosis model meihodology that 
predicts diameter growth for large trees and height growth for 
smaller trees. For the set of smaller trees, d.b.h. was not in- 
cluded in the list of possible predictor variables. The funda- 
mental reason for modeling the two size classes of trees dif- 
ferently is that for small trees, either diameter at breast height is 
not available, or is not as meaningful a dimension as it is for 
larger trees. 


Crown Width 


For 160 trees 3.5 inches (8.9 cm) d.b.h. and larger, the 
general model is: 


In(CW) = bo + b,yIn(D) + besln(H) + bgln(CL) 


In this and all subsequent models, a log-log model form was 
chosen in order to linearize the geometric relationships and 
stabilize the variance structure commonly encountered when 
working with a wide range of data (for instance, when variance 
is proportional to tree size). In this and subsequent models, the 
second assumption is borne out by the uniform random patterns 
shown in scatterplots of residuals versus predicted values in the 
log scale. 

Itis interesting to note that in this model no stand density term 
is included as a measure of competition. However, if the log 
model were inversely transformed to the natural scale, the 
combination of the CL term with its positive coefficient, and H 
with its negative coefficient (table 3) could be interpreted as 
crown ratio. The effect of competition on predicted crown width 
would be represented by crown ratio instead of stand density. It 
is assumed that this effect would be positive, that is, crown width 
would increase with increasing crown ratio. 

Table 4 is an analysis of variance showing the improvement 
in the regression sum of squares when variables in the general 
model are allowed to vary by species. For example, varying the 


Table 3.—Coefficients for estimating In (crown width) of trees 3.5 
inches d.b.h. and larger: 
In(CW) = bo, + b,In(D) + bzIn(H) + bzin(CL) 


Variable coefficients 


Intercept In(H) 
Species' bo, bo, 
BP —0.91984 — 0.07299 
WP 4.30800 — 1.37265 
WC 2.79784 — .89666 
WL 2.31359 — .80919 
WH 1.32772 — .52554 
AF 1.74558 = 392 
LP 1.06804 — .55987 
GF 2.20611 — .76936 
ES 3.76535 — 1.18257 
DF 3.02271 — 1.00486 
PP 1.62365 — .68098 
Variables Variable coefficients 
In(D) b; = 1.08137 
In(CL) bs = 0.29786 


33.37909/23 
6.66189/136 
0.04898 
0.8336 


Regression sum of squares/d.f. 
Error sum of squares/d.f. 

Mean square error 

Re 


‘See text page 2 for species list. 


intercept term by species improves the regression by a signifi- 
cant amount. Then, when either In(D) or In(H) are allowed to 
vary by species, the fit is further significantly improved. Of these 
two alternatives, the latter has a slight advantage. Although 
fitting In(CL) by species, or any more successively complex 
model also improves the regression, the signs of the coeffi- 
cients for the H and CL terms become inconsistent (some 
species with positive coefficients and some with negative coeffi- 
cients for a single parameter). The inconsistencies of sign 
among species could not be logically explained, for example, as 
a pattern indicating ranking of species by tolerance. Therefore it 
is assumed that the terms change sign because the model is 
overfitting the particular data set. Acceptable models were thus 
constrained to those with consistent coefficients. The final mod- 
el chosen is that in which the intercept and height terms vary by 
species (indicated by the subscript / on the parameter): 


In(CW) = bo, + biIn(D) + beIn(H) + b3in(CL) 


The coefficients for this model, shown in table 3, gave amean 
square error of 0.04898. The equation overpredicts crown width 
in the natural scale for four species and underpredicts for seven 
species, as shown in table 5. The mean crown width is under- 
estimated by 5.7 ft (1.7 m). 

Negative bias (mean residual deviation from zero), intro- 
duced when the inverse logarithmic transformation is used to 
convert log-normally distributed estimates to the original natural 
scale, can be approximately corrected by adding one-half the 
residual variance to the estimate on the log scale (Baskerville 
1972). This amounts to multiplying the estimate of crown width 
in the natural scale by exp[/2MSE]. Thus a factor of e 17-4899) = 
1.025 may be applied to the estimate in the natural scale to 
correct for underprediction. After “bias adjustment,” mean 
crown width is overestimated slightly, 1.2 ft (0.4 m). 


Table 4.—Analysis of improvement of In(crown width)' regression models attainable by varying parameters by species, trees 3.5 
inches d.b.h. and larger 


(1) 


(2) 


(4) 


(5) 


(6) 


(7) 


(8) 


Do 


Bo 


b,In(D) + baln(H) + bsln(CL) 


b,In(D) + baln(H) + bsln(CL) 


b,In(D) + bsIn(H) + b«ln(CL) 


—s 


b,In(D) + bsIn(H) + bgin(CL) 


b,In(D) + beln(H) + b3In(CL) 


b,In(D) + baIn(H) + bain(CL) 


b,In(D) + boln(H) + bgin(CL) 


b,In(D) + baIn(H) + b3in(CL) 


b,,In(D) + baIn(H) + bgIn(CL) 


Model? 


Source 


Remarks 


Degrees 
of 
freedom 


Marginal 


sum of Meansquare 


squares 


complete model 


Reduction due 
to model (1) 
Reduction due 
to model (2) 
-reduction due 
to model (1) 
Reduction due 
to model (3) 
-reduction due 
to model (2) 
Reduction due 
to model (4) 
-reduction due 
to model (2) 
Reduction due 
to model (5) 
-reduction due 
to model (2) 
Reduction due 
to model (6) 
-reduction due 
to model (4) 
Reduction due 
to model (7) 
-reduction due 
to model (3) 
Reduction due 
to model (8) 
-reduction due 
to model (4) 


Error 


‘Mean In(CW) = 2.53253 feet: standard deviation = 0.50183. 


2Terms varying by species are underlined. 
STabulated F values: 


Fi%6,001 = 2.49. 


inconsistent 
signs in CL terms 


inconsistent 
signs in H terms 


inconsistent 
signs in CL terms 


inconsistent signs 
in H and CL terms 


inconsistent signs 
in H and CL terms 


3 


116 


24.09107 


7.46691 


1.45889 


1.82111 


2.34274 


1.10706 


1.36227 


1.35748 


4.50798 


reduction 


8.03036 


.74669 


.14589 


18211 


.24327 


11071 


.13623 


.13575 


.03886 


F? 


206.65 


19.21 


6.03 


2.85 


3.51 


3.49 


Table 5.—Summary of residuals (natural scale) from prediction equations for 
crown width of trees. Bias adjustment factors are 1.025 for trees 
greater than or equal to 3.5 inches d.b.h., and 1.031 for trees less than 
3.5 inches 


Trees =3.5 inches d.b.h. Trees <3.5 inches d.b.h. 
> (obs-pred) > (obs-pred) 
No. of without with bias _No. of without with bias 
Species' trees correction correction trees correction correction 


BP 7 — SOFA 193 12 0.61 —0.19 
WP. 8 1.59 =) 3i5) 17 = 203 22 
WC 25 ke5 = 11:26 25 1.70 — 2.00 
WL 9 = 08} — 3.08 15 1.25 =e 
WH 12 15 — 5.66 16 2.08 20 
AF 10 — .80 — 3.00 15 58 =n 
He 7 5.68 ShTAS} 12 29 —7.08 
GF 22 1.80 —6.09 29 1.55 —2.40 
ES 8 2235) = 43 14 52 198 
DF 23 9.95 2.07 28 4.19 78 
PP 29 21.02 10.30 27 8.84 5.66 
Weighted mean 
residual (ft/tree) 5.67 = iss) 2.48 —0.07 
(m/tree) (1.73) (— 0.36) (0.76) (—0.02) 


‘See text page 2 for species list. 


Table 6.—Coefficients for estimating 
In(crown width) of trees less than 3.5 inches: 
In(CW) = b,, In(H) + bzIn(CL) + b3In(BA) 


The general model for prediction of crown width for 210 trees 

less than 3.5 inches (8.9 cm) d.b.h. is: 
Variable coefficients 
In(CW) = b,yIn(H) + beln(CL) + b3ln(BA) Species’ In(H) 
b 

In this model for smaller trees the density effect, as measured i 
by the basal area term, is positive rather than negative as might 
be expected (table 6). One possible explanation is that there is 


no competition effect among small trees due to stocking. In- BP 0.07049 
stead, the density effect is positive in that small trees in a WP 37031 
well-stocked stand expand laterally to utilize all available space. ye ue 
As the stand matures with time, closure of the canopy begins, WH 25622 
and individual tree crowns start to compete, stocking may nega- AF 33722 
tively affect crown width. This argument would assume that the LP 26342 
measured stands are mainly even-aged, and would not hold for GF 38503 
suppressed trees with weak crowns. That more than three- ES 33089 
fourths of the data are from dominant or codominant trees helps DF 32874 
support this explanation. PP .36380 
Variables Variable coefficients 
In(CL) bs = .28283 
In(BA) b, = .04032 
Regression sum of squares/d.f. = 371.02441/13 


Error sum of squares/d.f. 11.89093/197 
Mean square error 0.06036 
R? = 0.9689 


'See text page 2 for species list. 


Table 7 shows the effects on the general model of varying 
parameters by species. Since the intercept term in the original 
model was not significantly different from zero, all subsequent 
regressions were fitted through the origin. Varying In(H), In(CL), 
or In(BA) by species significantly improves the fit. However, 
inconsistency in signs for the In(CL) or In(BA) terms cccurs in all 
except model (2). Thus model (2) was chosen for predicting 
crown width for trees less than 3.5 inches (8.9 cm) d.b.h.: 


In(CW) = b,In(H) + beln(CL) + b3in(BA) 


where the subscript / denotes a different coefficient for the In(H) 
term for every species. 

The coefficients in table 6 gave a mean square error of 
0.06036. The model underpredicts crown width by an average 
of 2.5 ft (0.8 m) (table 5). A factor of e '7°°°) = 1.031 may be 
applied to the estimate in the natural scale to correct for nega- 
tive bias. When this is done, crown width is overpredicted by 
only 0.07 ft (0.02 m). 


Foliage Weight 


The general model for predicting needle weight developed 
from 114 trees 3.5 inches (8.9 cm) d.b.h. and larger is: 


In(WT) = Do + b,yIn(D) + bealn(H) + bgln(CL) + 
bgln(A) + bsln(TPA) + begln(DREL) 


The analysis of variance in table 8 is slightly different from that 
in tables 4 and 7. It was not possible to fit ine complete mode! 
(species x parameter interactions for seven variables) in the 
REX regression program because of the dimension of the prob- 
lem. Instead, the improvement in fit was analyzed between two 
models at a time by forming an F ratio in which the denominator 
is the mean square of the model with the most parameters in it. 
For instance, the general model was fit in (1). Then in (2) the 
intercept was varied by species, and the reduction in error 
between (1) and (2) associated with 9 d.f. was measured with an 
F-test. Since this reduction was significant, the species- 
intercept term was carried through in each successively more 
complex model, (3) through (8). Beyond fitting the intercept term 


Table 7.—Analysis of improvement of In(crown width)' regression models attainable by varying parameters by species, trees less 


than 3.5 inches d.b.h. 


Model? Source 


(1) b,In(H) + boIn(CL) + bsln(BA) Reduction due 
to model (1) 
Reduction due 
to model (2) 
-reduction due 
to model (1) 
Reduction due 
to model (3) 
-reduction due 
to model (1) 
Reduction due 
to model! (4) 
-reduction due 
to model (1) 
Reduction due 
to mode! (5) 
-reduction due 
to model (2) 
Reduction due 
to model (6) 
-reduction due 
to model (2) 
Reduction due 
to model (7) 
-reduction due 
to model (3) 


(2) by In(H) + beln(CL) + bsln(BA) 


(3) byIn(H) + bsIn(CL) + bsln(BA) 


(4) b,In(H) + baIn(CL) + bgIn(BA) 


(5) by In(H) + bzIn(CL) + b3ln(BA) 


(6)  b,\In(H) + baln(CL) + bIn(BA) 


— 


(7) b;In(H) + beIn(CL) + bgIn(BA) 


complete model 
by In(H) + bzIn(CL) + b3In(BA) Error 


‘Mean In(CW) = 1.26648 feet; standard deviation = 0.46956. 
Terms varying by species are underlined. 
STabulated F values: F}9;995 = 1.63. 


Degrees Marginal 
of sum of Mean square 


Remarks freedom squares reduction F? 


3 364.13574 121.37858 2063. 


10 6.88867 .68887 EVAL 
Inconsistent 
signs in CL terms 10 6.89868 .68987 tiewS 
Inconsistent 
signs in BA terms 10 6.77222 .67722 11.51 
Inconsistent 
signs in CL terms 10 1.16358 .11636 1.98 
Inconsistent 
signs in BA terms 10 .81690 .08169 1.39 
Inconsistent 
signs in BA terms 10 1.07251 10725 1.82 


Inconsistent signs 
in CL and BA terms Were 


10.41343 .05883 


Table 8.—Analysis of improvement of In(foliage weight)' regression models attainable by varying parameters by species, trees 


3.5 inches d.b.h. and larger 


Model? Source 


(1) bo + byln(D) + beln(H) + bgln(CL) Reduction due 
to model (1) 
Residual from 
model (1) 
Reduction due 
to model (2)-(1) 
Residual from 
model (2) 
Reduction due 
to models (3)-(2) 
Residual from 
model (3) 
Reduction 

to models (4)-(2) 
Residual from 
model (4) 
Reduction due 
to model (5)-(2) 
Residual from 
model (5) 
Reduction due 
to models (6)-(2) 
Residual due 

to model (6) 
Reduction due 
to models (7)-(2) 
Residual from 
model (7) 
Reduction due 
to models (8)-(2) 
Residual from 
model (8) 


+ baln(A) + bsln(TPA) + beln(DREL) 
(2) bo, + b,In(D) + baln(H) + bgln(CL) 

+ baln(A) + bsln(TPA) + beln(DREL) 
(3) Do, - b,In(D) + beln(H) + bgln(CL) 

+ b4ln(A) + bsln(TPA) + beln(DREL) 
(4) Do + b,In(D) + b2In(H) + bgln(CL) 

+ byln(A) + bsln(TPA) + b,eln(DEL) 
(5) Do + b,In(D) + boln(H) + bg In(CL) 

+ bgln(A) + bsln(TPA) + beln(DREL) 
(6) bo + byIn(D) + beln(H) + baln(CL) 

+ b4aln(A) + bsIn(TPA) + beln(DREL) 
(7) bo + byIn(D) + baln(H) + bgln(CL) 

+ byln(A) + bsin( TPA) + Degin(DREL) 
(8) bo + byIn(D) + beln(H) + baln(CL) 


+ baln(A) + bsIn(TPA) + bein(DREL) 


'Mean in(W7) = 3.41782 Ib; standard deviation = 1.24050. 
?Terms varying by species are underlined. 
3Tabulated F values: F830, = 2.59; F49005 — 198 


by species, there was no significant species interaction with any 
of the variables in the general model except for In(CL). How- 
ever, since the signs of the coefficients for the species x In(CL) 
term were inconsistent, this model was rejected as possibly 
overtitting the data. 

Thus, the final model is that in which just the intercept term 
varies by species: 


In(WT) = bo, + biIn(D) + beln(H) + b3ln(CL) + 
bgln(A) + bsIn(TPA) + beln(DREL) 


Remarks 


Degrees Marginal 
of sum of 
freedom squares Meansquare F? 


6 152.41829 25.40305 126.59 


107 21.47126 .20067 
9 8.43661 .93740 7.05 
98 13.03464 .13301 
9 1.35027 .15003 1.14 
89 11.68435 .13128 
9 1.60567 .17841 1.39 
89 11.42897 12842 
9 2.20900 24544 2.02 
inconsistent 
signs in CL terms 89 10.82565 12164 
fe) .75875 .08431 61 
inconsistent 
signs in A term 89 12.27590 .13793 


9 81401 .09446 .66 


inconsistent 
signs in TPA term 89 


12.22063 13731 


9 1.41581 .15731 1.20 


inconsistent 
signs in DREL terms 89 


11.61885 * .13055 


The coefficients in table 9 gave a residual mean square of 
0.13301. Note that since there were no data available for lodge- 
pole pine, no parameters were estimated for this species. Table 
10 shows for which species the model overpredicts and under- 
predicts. The mean prediction underestimates foliage biomass 
by 86.9 lb (39.4 kg). Acorrection for bias of e'*'*°" = 1.069 may 
be applied, resulting in an overprediction of 5.7 |b (2.6 kg) per 
tree after adjustment. 

In the above model, predicted foliage weight decreases as 
trees per acre increase, an indication of competitive stress due 
to stocking. The positive crown length and negative height, 
which reflect crown ratio in the natural scale, and positive rela- 
tive diameter terms may all be interpreted as measures of an 
individual tree’s competitive status. In addition, the negative 
age term coupled with the positive diameter term may be inter- 
preted as an estimator of mean annual diameter increment in 
the natural scale, which would be expected to vary directly with 
foliage weight. As tree age is not always a readily available 
measurement in inventory data, an alternative equation for 
foliage weight involving a transformation of the diameter and 
age terms is presented in a later section. 


Table 9.—Coefficients for estimating In(foliage weight) 
of trees 3.5 inches and larger: 
In(WT) = bo, + b,In(D) + bzin(H) + b3in(CL) 
+ b,In(A) + bsIn(TPA) + beln(DREL) 


Variable coefficients 


In developing an equation for foliage weight prediction for 210 
trees less than 3.5 inches (8.9 cm) d.b.h., the general model 
form is: 


Intercept 
Species’ bo, 
BP 2.62251 
WP 2.66607 
WC 3.05935 
WL 1.75654 
WH 2.65457 
AF 3.06017 
GF 3.11508 
ES 3.30085 
DF 2.70587 
PP 2.45249 
Variables Variable coefficients 
In(D) b, = 2.08624 
In(H) bs = —1.07705 
In(CL) b3 = .69082 
In(A) b, = — -30885 
In(TPA) b; = — -14210 
In(DREL) on = 39924 
Regression sum of squares/d.f. = 160.85490/15 
Error sum of squares/d.f. = 13.03465/98 
Mean square error = (aale'stoh 
R? ~ 0.9250 


‘See text page 2 for species list. 


In(WT) + bo + byIn(CL) + bein(A) + bgln( TPA) 


The analysis of variance of species by parameter interactions 
for this model is shown in table 11. Varying the intercept and 
In(CL) by species resulted in significant improvement while 
maintaining consistency in signs between species. Improve- 
ments due to varying other parameters by species were also 
significant, but integrity of sign was no longer maintained. 
Therefore, the model chosen for predicting foliage weight of 
trees less than 3.5 inches (8.9 cm) is: 


In(WT) = bo + byIn(CL) + boin(A) + bsln( TPA) 


where the subscript / denotes a parameter varying by species. 


Table 10.—Summary of residuals (natural scale) from prediction equations for 
foliage weight of trees. Bias adjustment factors are 1.069 for trees 
greater than or equal to 3.5 inches, and 1.126 for trees less than 3.5 


inches 


Trees =3.5 inches d.b.h. Trees <3.5 inches d.b.h. 
> (obs-pred) > (obs-pred) 


No. of without 


with bias No. of without with bias 


Species' trees correction correction trees correction correction 


BP = 9.44 
WP 2 SOIT 
WC 22 — 218.20 
WL 3 — 05 
WH o toi 
AF 10 — 2:22 
EP 0) — 

GF 17 48.58 
ES 3 4.16 
DF 20 — 52.94 
PP 29 515.50 


Weighted mean 
residual (Ibs/tree) 86.89 
(kg/tree) (39.41) 


=3:92 12 1.31 — 0.07 

— 4.65 17 1.85 11-90 

— 330.13 25 5.77 55 
= .98 15 —16 — 1.90 
=5:23 16 1.34 = 2:96 

— 30.10 15 a0 .28 
— 12 0.01 — 1.26 
Eee 29 3.88 —§.35 

— 5.42 14 1.56 —3.13 
—172.73 28 7.61 1.43 
362.32 27 Sal 1.10 
—5.69 3.71 al an 
(— 2.58) (1.68) (—0.51) 


‘See textpage 2 for species list. 


Table 11.—Analysis of improvement of In(foliage weight)' regression models attainable by varying parameters by species, trees 


less than 3.5 inches d.b.h. 


Degrees Marginal 
of sum of Meansquare 
Model? Source Remarks freedom squares reduction F? 

(1) Do + byIn(CL) + bealn(A) + bsln( TPA) Reduction due 

to model (1) 3 259.86645 86.62215 413.06 
(2) Do, + b,In(CL) + bsln(A) + bs3ln( TPA) Reduction due 

to model (2) 

-reduction due 

to model (1) 10 32.33326 3.23326 15.42 
(3) Do + b,In(CL) + bein(A) + bsln( TPA) Reduction due 

to model (3) 

-reduction due 

to model (2) 10 5.37841 53784 2.56 
(4) Do, + b,In(CL) + baIn(A) + bsln( TPA) Reduction due 

to model (4) 

-reduction due inconsistent 

to model (2) signs in A term 10 6.87378 .68738 3.28 
(5) bo + b,In(CL) + bsln(A) + bgln(TPA) Reduction due 

to model (5) 

-reduction due inconsistent 

to model (2) signs in TPA term 10 5.20996 .52100 2.48 
(6) Do, o b,In(CL) - ba|n(A) + b3ln( TPA) Reduction due 

to model (6) 

-reduction due inconsistent 

to model (3) signs in A term 10 5.54004 .55400 2.64 
(7) Do + b,In(CL) + baln(A) + bgln( TPA) Reduction due 

to model (7) 

-reduction due inconsistent signs 

to model (3) in TPA terms 10 5.38501 .53850 PAN 
(8) Do, + b,In(CL) + baIn(A) + b3In( TPA) Reduction due 

to model (8) 

-reduction due __ inconsistent signs 

to mode! (4) inAand TPAterms 10 3.79663 .37966 1.81 

complete model 
Do + b,In(CL) + ba|n(A) - bs In( TPA) Error inconsistent signs 
inAand TPAterms 166 34.81201 .20971 


‘Mean In(WT) = —0.075104 Ib; standard deviation = 1.27874. 
“Terms varying by species are underlined. 
STabulated F values: Fi8.o05 = 1.89; Fi26001 = 2.43. 


Mean square error for the above model is 0.23751, and 
coefficient values are shown in table 12. Table 10 summarizes 
the residuals for the model. Foliage weight is underestimated by 
an average of 3.7 Ib (1.7 kg) per tree. Negative bias may be 
compensated for by multiplying the estimate in the natural scale 
by e'7*75) = 1.126, giving an overprediction of 1.1 Ib (0.5 kg) 
per tree after adjustment. 

The present approach for estimating foliage biomass uses 
two model forms; one for small trees and one for large trees — 
each of which is applicable to all species. The data used in 
developing these equations were a subset of Brown's (1978) 
data. His approach was to fit one or more separate model forms 
for total crown biomass (foliage and branchwood) for trees 
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greater than 1 inch d.b.h. and for trees less than 2 inches d.b.h., 
by species and crown class. He further fitted separate regres- 
sion equations for predicting the proportion of total crown 
biomass in foliage and branchwood components. Foliage 
weight can be estimated indirectly using Brown’s method by 
applying these two equations in succession to the data. The 
following method was used to compare the present approach 
with Brown’s approach for predicting foliage biomass. The best 
combination (that yielding lowest mean square error) of Brown's 
crown biomass and foliage proportion equations for a species 
was solved for each tree in the data set used to develop the 
current equations. Residual sums of squares [{(observed — 
predicted needle weight)*] were then compared for the two 
methods. Results are as follows: 


Table 12.—Coefficients for estimating In (foliage weight) of trees 
less than 3.5 inches d.b.h.: 
In(WT) = bo, + b,In(CL) + bzIn(A) + b3In(TPA) 


Variable coefficients 


Intercept In(CL) 
Species’ bo, b,, 
BP —A.Gioile 1.47513 
WP — 2.15894 1.48969 
WC — 2.64034 1.69973 
WL — 5.02156 2.31835 
WH — 4.22701 2.22534 
AF — 2.03919 1.64942 
EP — 3.38394 1.96060 
GF — 2.78090 1.90272 
ES — 3.30673 2.27613 
DF — 2.30430 1.52896 
Pe — 3.02050 1.88712 
Variables Variable coefficients 
In(A) bs = .22823 
In( TPA) b, = — .13550 


Regression sum of squares/d.f 
Error sum of squares/d.f. 
Mean square error 

R2 


= 297.57812/23 

44.17596/186 
0.23751 
0.8707 


Il 


Il 


"See text page 2 for species list. 


Residual sums of squares for prediction 
of In (foliage weight) 


Trees <3.5 Trees 23.5 
inches d.b.h. inches d.b.h. Total 
Current 44.17596 13.0346 57.21056 
models 
Brown's 69.9460 10.08098 80.0269 
models 


The current model for trees greater than or equal to 3.5 inches 
(8.9 cm) d.b.h. had a larger residual sum of squares than that 
obtained from Brown’s methods. For trees less than 3.5 inches 
(8.9 cm), residual sum of squares obtained by the present 
method are lower. Overall there appears to be no loss in predic- 
tive ability by using a single model form as presented here. 
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Alternative Models for Foliage 
Weight Prediction 


The prediction models for both large and small trees use the 
variable In(age), which is not available in the prognosis system 
(Stage 1973) in which the equations are to be used. Thus the 
following alternative models are presented. 

For trees 3.5 inches (8.9 cm) and larger, the alternative model 
is a transformation of the general model involving the diameter 
and age terms. Periodic (diameter)* increment, which is avail- 
able in the prognosis model, is substituted for mean annual 
(diameter)* increment in the following transformation: 


In(WT) = Do + b,In(D) + baln(H) + bgln(CL) + 
baln(A) + bsin(TPA) + bgln(DREL) 
(general model in log linear form) 


or, 
WT = e%)- D1 - HP2 - (CL)P3 - AP4 - (TPA)®5 - (DREL) 
(general model in exponential form) 
Multiplying by D*°4 
D??4 
(b, + 2b ) 
4 


wr=e%-D SH 


(CL)°3 - (TPA)°S - Ene 
which, in logarithmic form, is 
In(WT) = bo + (by + 2b,)-In(D) — baln (DO 
+ beln(H) + b3ln(CL) + bsln(TPA) + SNORER 
and substituting periodic [A (D*)] for mean annual [D*/A)] incre- 
ment, 


In(WT) = bo + (b; + 2b.)In(D) — b,ln[A(D?)] + 
boIn(H) +b3In(CL) +bsIn(TPA) + bsln(DREL) 
(alternative model) 


Kittredge (1944) noted that it would be expected that “the 
amount of foliage could be estimated from the periodic annual 
growth because the increment of stem wood is determined by 
the amount of foliage which is carrying on photosynthesis in that 
period” (p. 906). It is important to note that the alternative model 
has not been reparameterized with the periodic annual incre- 
ment (PAI) as an independent variable since PAI was not mea- 
sured. It is merely a transformation of the original parameters 
and should be presented with a statement of caution about the 
relationship between the mean annual increment (MAI) and 
PAI. 

The substitution of periodic (diameter)* increment for mean 
annual increment in the alternative equation assumes that the 
relationship between (diameter)* and age is linear over a wide 
range of ages. The assumption breaks down when PAI and MAI 
diverge as typically happens over time except at the culmination 
of MAI. Figure 1 shows the difference in predicted foliage weight 
between the current and alternative models that can be ex- 
pected when PAI departs from MAI. When PAI is 100 percent 
greater than MAI, predicted foliage biomass increases 24 per- 
cent. PAI 60 percent of MAI produces a 25 percent negative 
change in predicted foliage biomass. 
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Figure 1—Change in predicted foliage biomass as periodic annual increment departs from mean annual 


increment. 


To compare whether the alternative model approximates the 
fit obtained by the general model, both models were applied to 
an independent stand data set* in which both periodic and mean 
annual increments were measured. The tree ages in this stand 
ranged from 2 to 560 years. For 95 trees in the test data set, the 
chi-square value for the sum of squared differences between 
needle weight prediction by the alternative model and prediction 
by the general model was 4.10 compared to the tabulated value 
of x? 95. 0005 = 134.24. The alternative model overpredicts 
needle weight outside of the range of data used to develop the 
general model (the sum of differences between predictions by 
the alternative and current models was 741.3 Ib [336.2 kg] for 
the 95 trees). However, with the low chi-square value, the 
alternative model may be judged adequate, and also uses 
independent variables that are more readily available. 

For trees less than 3.5 inches (8.9 cm) d.b.h., the following 
reparameterized alternative model substituting height for age is 
presented: 


In(WT) = bo, + byIn(H) + beIn(CL) + bgln(TPA) 


Mean square error is 0.24759 for the coefficients given in table 
13. This model underestimates needle weight for 9 of 11 spe- 
cies (table 14). On the average, it underpredicts needle weight 
by 3.6 Ib (1.6 kg). A correction for underestimation of e'*\*""" = 
1.132 may be applied giving a mean overprediction of 1.0 lb (0.5 
kg) after bias adjustment. Residual sum of squares obtained for 
this alternative equation is 46.0511, which Is also lower than 
that obtained from Brown’s method. Thus, the alternative 
formulation is judged acceptable. 


*St. Joe National Forest, Idaho. 
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Table 13.—Coefficients for estimating In (foliage weight) of trees 
less than 3.5 inches d.b.h., alternative model: 
In(WT) = bo, + biIn(H) + bzIn(CL) + bgln(TPA) 


Variable coefficients 


Intercept In(CL) 
Species’ Bo, bo, 
ee 
BP — 2.63387 1.35092 
WP — 1.94951 1.22023 
WC — 2.24876 1.37600 
WL — 4.73762 1.98479 
WH — 4.17456 2.00749 
AF — 1.60998 1.32649 
EP. — 3.13488 1.62368 
GF — 2.43200 1.60270 
ES — 2.93508 1.96125 
DF — 2.05828 1.25837 
PP — 2.74410 1.58171 
Variables Variable coefficients 
In(H) b, = .40350 
In( TPA) bs = —.112975 


295.70288/23 
46.05106/186 
0.24759 
0.8653 


Regression sum of squares/d.f. 
Error sum of squares/d/f. = 
Mean square error 
R? 


| 


| 


'See text page 2 for species list. 


Table 14.—Summary of residuals (natural scale) from alternative 
prediction equation for foliage weight of trees less than 
3.5 inches d.b.h. The bias correction factor is 1.132 


___% (obs ~ pred) __ 


No. of without with bias 
Species’ trees bias correction correction 
BP 12 1.41 0.11 
WP vA 1.85 — 1.74 
WC 25 5.74 74 
WL 15 —.45 —215 
WH 16 1.14 —3.00 
AF 5 5.83 .87 
LP 12 —.09 Soe 
GF 29 3.44 — 5.29 
ES 14 1.63 =2:85 
DF 28 7.39 1.46 
PP 27 5.23 1.36 
Weighted mean residual 
(Ib/tree) 3.62 =1.03 
(kg/tree) (1.64) (—.47) 


‘See text page 2 for species list. 
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Relationships between tree crown width and diameter at breast height, height, 
crown length, and basal area per acre, and between tree foliage biomass and 
diameter at breast height, height, crown length, age, relative diameter, and trees per 
acre are presented for 11 conifer species in the Northern Rocky Mountains. 


KEYWORDS: crown width, foliage biomass 
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The Intermountain Station, headquartered in Ogden, 
Utah, is one of eight regional experiment stations charged 
with providing scientific knowledge to help resource 
managers meet human needs and protect forest and range 
ecosystems. 

The Intermountain Station includes the States of 
Montana, Idaho, Utah, Nevada, and western Wyoming. 
About 273 million acres, or 85 percent, of the land area in the 
Station territory are classified as forest and rangeland. These 
lands include grasslands, deserts, shrublands, alpine areas, 
and well-stocked forests. They supply fiber for forest in- 
dustries; minerals for energy and industrial development; and 
water for domestic and industrial consumption. They also 
provide recreation opportunities for millions of visitors each 
year. 

Field programs and research work units of the Station 
are maintained in: 


Boise, Idaho 


Bozeman, Montana (in cooperation with Montana 
State University) 


Logan, Utah (in cooperation with Utah State 
University) 

Missoula, Montana (in cooperation with the 
University of Montana) 


Moscow, Idaho (in cooperation with the Univer- 
sity of Idaho) 


Provo, Utah (in cooperation with Brigham Young 
University) 


Reno, Nevada (in cooperation with the University 
of Nevada) 
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